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Abstract 

We establish the equivalence between the quantum evolution of spatially 
homogeneous oscillations of a scalar field and that of an analogous classical 
system with certain random initial condition. We argue that this observation 
can be used for numerical simulation of the Universe in the preheating epoch. 
We also explicitly demonstrate that the phenomenon of parametric resonance 
that leads to preheating is simultaneously an effective mechanism for generating 
quantum decoherence of the Universe. 
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1. Calculation of the reheating temperature is one the most important issue in 
inflationary cosmology. The temperature achieved soon after the inflation epoch is 
the maximal temperature in the whole history of the Universe and must be sufficiently 
large for the generation of baryon number to be effective. According to the inflationary 
picture, the Universe is reheated by the process of transferring energy out of spatially 
homogeneous oscillations of the scalar infiaton field. The old theory of the reheating 

describes this energy transfer as occurring through the decay of the infiaton quanta 
and predicts a relatively low speed for this process in theories with small couple 
constants. However, a recent revision shows that in many case there exists a more 
rapid channel of energy transfer, which is closely related to the existence of parametric 
resonance in the spectra of the fields interacting with the infiaton and of the infiaton 
itself. These parametric resonance modes, typically forming continuous energetic 
bands, grow exponentially on the background of the infiaton field oscillations, and at 
some moment begin to take a large amount of energy from these oscillations. It has 
been argued in ref. that due to the exponential character of the growth, the decay 
of the infiaton field oscillations is a rapid process of "explosion" rather than a slow 
damping. The amplitude of these oscillations drops to a small value after some time 
interval of the same order as of magnitude of the period of oscillations. This new 
mechanism, named in ref. "preheating" , may change considerably the prediction of 
the reheating temperature. It has been argued that the phenomenon of parametric 
resonance may be important in other aspects as well 0, ^. 

The complete theoretical description of preheating is still lacking right now, due 
to the complexity of interaction between the modes that appear via parametric reso- 
nance and the background, as well as between the modes themselves. However, some 
important steps have been made in this direction (see, for example 0, 0). In one of 
the approaches ^ the problem is treated in the Hartree-Fock approximation, where 
one takes into account the back-reaction of the parametric-resonance modes to the 
background but neglects the scattering between these modes. The result is quite 
surprising: while it confirms that parametric resonance is the main mechanism of the 
damping of the infiaton field in the first time period, in the long run the amplitude 
of the oscillations typically remains constant (in some cases it is not much smaller 
than the initial amplitude before the explosion), or at least decays very slowly. This 
shows that it takes much longer time, of many order of magnitudes larger than the 
period of oscillations, to make the amplitude goes to 0. At late times, the created 
particles are mostly of very small energy, which means that the thermalization pro- 
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cess may also take a longer time interval than has been expected. The calculations 
in ref. ^ , however, rely heavily on the Hartree-Fock approximation, so there remains 
the question whether these results reflect the real situation. 

In this paper we do not attempt to build a complete theory of preheating. How- 
ever, we will point out a possible way to perform reliable calculations in the nonlinear 
regime after inflation. Namely, we show that the behavior of the quantum system 
under consideration, during the preheating epoch, is equivalent to the evolution of 
the classical field, when one randomizes over a particular set of initial conditions. 
This result is established by comparing the perturbative series of the two theories. 
This fact makes possible numerical simulations of the Universe after inflation, since 
the problem is now completely classical. As a by-product, we demonstrate explicitly 
the development of decoherence of the Universe. We show that the phenomenon of 
parametric resonance provides an effective mechanism for the latter. 

2. For simplicity we will work in the model of one scalar field only. In other words, 
we neglect the interaction of the infiaton with other particles and take into account 
its self-interaction only. The Lagrangian of the model is taken in the standard form. 

We will also neglect the expansion of the Universe and work in the ordinary Minkowskian 
space-time. At the moment when inflation completes, the infiaton field has a nonzero 
expectation value = 0o which is a constant over the whole space. Soon after that 
the field begins oscillations around = 0. For convenience the following picture is 
used in further discussions: at t < there is an extra force term Jcf) in the Lagrangian 
so that the expectation value of is a nonzero constant 0O; but at t = this source 
is suddenly turned off and the field starts oscillate. We will be interested in the time 
dependence of the mean value of 0, which will be denoted as 0(t). 

The linearized classical equation for the modes of momentum k of the field, (f){t, k), 

dl + ujl + ^<Pl(t))<Pit.^) = ^ (1) 

has the same form as that of the oscillator whose frequency is a periodic function of 
time, uj"^ = uj1 + It is well known that at some values of uj^ this system exhibits 

parametric resonance: there does not exist any solution to eq.(0) which remains 
finite on the whole time axis |]^. A typical solution to eq.(|l|) grows exponentially 
in the limits t —>■ +oo and/or t —oo. This phenomenon occurs when ujk lies in 
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certain energetic bands ("resonance bands"). We will consider the situation of broad 
resonance, which correspond to large values of 0o, 4>o ~ or larger. 

In the quantum case, the existence of parametric resonance leads to the amplifi- 
cation of the quantum fluctuations of the modes with uJk inside the resonance bands. 
When the amplitude of these fluctuations are still small, one can make use of the 
standard Bogoliubov transformations to find the quantum state of the system [Q. 
However, the interesting regime is that of late times where the non-linearity becomes 
essential, and one needs another approach to attack the problem. 

The general method to find real-time evolution of quantum fields is the Schwinger- 
Keldysh closed-time-path formalism p, |^, |10[ . At diagrammatic level, the formalism 



gives rise to the same set of Feynman rules as that for calculating the S'-matrix, 
except that every internal integral in Feynman diagrams is now performed along a 
contour that goes from t = — oo to t = +00 and then goes back to t = —00 (Fig.l), 
and instead of one Feynman propagators G{x, y) there exists four ones depending on 
whether xq and yo lie on the upper or the lower parts of the contours, 

G++(a;,y) = (O|T0(x)0(y)|O), G+_(x,y) = (O|0(y)0(x)|O), 

= {0\<Pix)cj>{y)\Q), G^4x,y) = (0|f 0(a;)0(|/)|O), 

It is easy to see that these propagators arise from the propagators on the contour C, 
Gc{x,y) = {0\Tc4>{x)(f){y)\0) , where Tc is the notation for time-ordering along C. 
Note also that the four propagators are equal at x = y. 

To derive the Feynman rules, we will decompose the quantum field into a classical 
part (poit) and quantum fluctuation (p, 

4>it) = Mt) + kt) 

where (f)o{t) = (pQ at t < and at t > satisfies the field equation 



dl(j)o + ml(j)o + X(f)l = 



The Lagrangian for is then 



and corresponds to the Feynman rules shown in Fig. 2. 

Before turning to the calculation of diagrams, let us consider in more details 
the propagators on the background (po(t)- Consider, for instance, G++{x,y). It is 
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convenient to use the mixed representation (coordinate in time and momentum in 
space) where this propagators can be written into the form, 

r r/k 

G++(x,i/) = y^e^'^(''-y)G++(a;o,2/o;k) 

where 

iG++(xo, z/o; k) = ^(xo - l/o)/f (a;o)/2 (l/o) + ^(l/o - a;o)/2 (a;o)/i^i/o) (2) 

We have introduced the mode functions 2 which are the two hnearly independent 
solutions to the equation 

dl + ^l + \<Pl(t))fUt)=^ (3) 
with the following boundary conditions at t < 0, 

where Uk = \/ + rn^ + The mode functions are normalized so that f[{t)f2{t) — 
fi{t)f^{t) = -i. If (j)o{t) is 0, /f and /| would be equal to e'^-^^Vv^ and e''^>=^ 
respectively. It is easy to see that = ifiit))*. 
The formulas for other propagators are similar, 

iG+_(xo,?/o;k) = f2{xo)fi{yo) 

iG_+(xo,yo;k) = fi{xo)f^{yo) 

iG^4xo,yo; k) = e{xo - yo)f2Mfi{yo) + ^(2/0 - a;o)/f(a;o)/2 (2/0) (4) 

Consider now a value of k where there is parametric resonance in eq.(|^). When 
parametric resonance is present, eq.(^) possesses two real solutions and f^{t) 

which satisfy the following conditions at t > 

/^(t + T) = e^^^/^(t) 

fy + T) =e-''=^/!(t) 

where is some real, positive constant depending on uj^- In particular, grows as 
t —>■ +00 while the behavior of /_ is opposite. The previously defined mode functions 
/i and /2 are linear combinations of /+ and /_. Not interested in the particular form 
of the coefficients, we write, 

A'{t) = auflit) + f,(t) = alfld) + Plftit), (5) 
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where we have made use of the fact that /i and /2 are complex conjugate and /+ 
and /_ are reaL Since as /+ is growing with t while /_ is falling, at large t the part 
f+(t) dominates over f-(t). If in eq.(|^) one neglects /_ and leaves only /+, there 
is no difference between /i and /2, except for an overall factor. If one also neglects 
modes outside resonance bands, which is natural since these modes are not enhanced, 
and substituting eq.(^ to eqs.(0), one finds that the four propagators are equal 
to each other to the leading order, 

iG++{x,y) = iG+^{x,y) = iG^+{x,y) = iG^^{x,y) = 

r r/k 

= ^G\x,y)= —-\a,\'flixo)fJyo) (6) 



. (27r)3 

Now let us consider the one-loop contribution to (p. The only diagram that con- 
tributes to is the tadpole one shown in Fig. 3a. This diagram is equal to 

<t^ii.^) = y J dyoj dyGc{x,y)(l)oiyo)G{y,y) = — J dy GR{x,y)(j)o{yo)G{y,y) (7) 
c 

where we have dropped the indices of G{y,y) since all the four Green functions are 
equal at coinciding points and introduced the retarded Green function, Gr = 6*++ — 
One can expect that at late Xq the integral in eq.(|^) is dominated by large 
values of yo where G{y,y) can be replaced by the leading term, eq.(^. So one finds, 

M^o) = J dyoGR{xo,yo;0)Myo) J |«fc/^(l/o)P (8) 

Eq.(^ can be represented graphically as in Fig.Sb, where each external line ending 
with a bullet is associated with the factor |afc|/+(?/o) if k is a resonance mode and in 
the opposite case. Note that grows exponentially with t, so 0i is also growing. At 
t when 01 becomes comparable to 0o oiie can expect that all terms of the perturbative 
series are of the same order, and the perturbation theory breaks down. We will try 
to extract the main contribution from each order of the perturbation theory in this 
regime. We will not be interested at latter times when 0i ^ 0o- 

Let us move to two-loop diagrams. Consider, for example, the one depicted in 
Fig.4a, 

iX^ r dk 

M^o) = -Y J dyodzo j^-^ Gcixo,yo]Q)(j)o{yo)Gc{yo, zo;k)Gciyo, zo;k)(f)o{zo)(l)iizo) 

(9) 
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If one is interested large values of xq, it can be expected that the important region of 
integration in r.h.s. of eq.(|) is that of large yo and zq. So, one may suspect that the 
integral can be calculated by replacing Gc{xo, yo; k) by its leading order contribution, 
eq.(||). One obtains after that 

M^o) = — — J dyodzo——^Gcixo,yo;0)Myo)if^iyo)yifU^o)fMzo)Mzo) 
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However, it is easy to see that the integration over zo gives zero, since it goes along 
the contour C and the integrand, (/^(zo))^0o(^o)0i(^o)5 is the same on the upper 
and lower parts of the contour. So, one should turn to the next-to-leading order. 
Denoting the next correction to eq.(|D as G^{x, y) = Gc{x, y) — G°(x, y), one writes, 

f dk 

02(a;o) = dyo dzo j^-^ Gc{xo, yo; O)0o(l/o)G'e(|/o, zq; k)/+(|/o)/+(?/o)0o(^o)0i(^o) 
r dk 

= X^J dyodzoj^-^Gc{xo,yo]0)(j)o{yo)Gc{yo,zo;k)fl{yo)fl{yo)Mzo)Mzo) (10) 

where in the last equation we have replaced G^j by Gc, making use of the fact that the 
leading term vanishes. Eq. ([I0|) can be rewritten in a simpler form using the retarded 
propagator Gr = G++ — 6*+-, 
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M^o) = J (^Vo dzo Gr{xo, yo; O)0o(l/o)Gij(?/o, zo; k)/+(yo)/+(l/o)0o(2o)0i(^o) 

— oo 

where the integrations over dyo and dzo are performed from — oo to +00. Eq.(|l^) can 
be represented graphically as in Fig.4b. As in the case of the one-loop graph, the 
calculation of the loop diagram, thus, is reduced to evaluation of a tree connected 
graph. 

This technique can be generalized to deal with other diagrams and higher orders of 
perturbation theory. Summarily, the technique is the following. To calculate a given 
Feynman diagram with I loops, one cuts / internal propagators, so that the obtained 
diagram is tree, and connected. If there are many ways to perform this cutting, one 
should take a sum over all possibilities. Each propagator that has been cut is associ- 
ated with two factors of attached to the two vertices that the former propagator 
had been connecting. The propagators that have not been cut are associated with 
the retarded propagator Gr. The obtained diagrams are computed by the standard 
way, where the integral are performed along the usual time axis. We note that this 
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reduction of loop graphs to tree ones is very similar to that found in another setting 
p!T| . We emphasize that our method is applicable in the regime when all terms of the 
perturbation series are of the same order. At latter times when 0i ^ 0o our approach 
does not work, since the pieces that have been neglected may become important. 

3. The fact that the diagrams obtained after the cutting procedure are tree 
points to the possibility of classical description of the problem. To find the latter, let 
us consider the Cauchy initial problem with the classical field equation 

{dl + m^)(f) + X(f)'' = (11) 

with the following initial condition at the asymptotics t — oo, 

0(t,x) = 0o(t) + /(^^^/^(^K'^ (12) 

where Ck are some arbitrary set of complex numbers. For to be real, we require that 
Ck = C-k- Since at t ^ —oo f^(t) is small and satisfies eq.(^, eq.([T^) is the solution 
to the field equation up to the first order of /+. The corrections to eq.(0), which start 
at the second order of /+, can be calculated iteratively from the field equation and 
the result can be represented in the form of tree Feynman diagrams. For example, 
the second-order contribution is represented by the first diagram in Fig. 5 together 
some third- and fourth-order graphs. In these graphs each external leg ending on a 
bullet is associated with the factor Ck/+, and each internal line corresponds to the 
retarded propagator. 

Despite the fact that these diagrams are similar to that obtained by the cutting 
procedure described above, there are two main differences between the two cases. 
First, the diagrams obtained by cutting loop graphs always contains an even number 
of external bullet legs. Second, the external bullet legs in the first case can be grouped 
into pairs so that the momenta running along legs of the same pair are of opposite 
signs, while in the graphs coming from eq.(|l2D the only requirement is that the sum 
of the momenta of the bullet legs is 0. Nevertheless, it is easy to show that there 
is a simple relation between the two cases. Namely, if one denotes as 0[i(:,x, Ck] the 
solution to the field equation with the boundary condition (|12|), then the quantum 
average of (p in our problem is given by an integral over Ck with a Gaussian measure, 

0M=/l>c,exp(-/^i^)#.x.cJ (13) 

(the integral in the r.h.s. of eq.(|13|) does not depend on x). To see this, one notices 
that the integration over Vck leaves only diagrams with even number of bullet legs 
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which can be grouped in paired of opposite momenta. The integration also adds a 
factor of to each pair with momentum (k, — k). So, the graphical representation 
of the r.h.s. of eq.(]TB|) is the same as those obtained by cutting the loop graphs. 

In eq.(|T3|), Ck are supposed to be real, Gaussian-distributed around with the 
standard deviation \ak\ and the Ck with different k are not correlated. Eg. ([131) is not 
the only possible, one could write, for example 

0(t) = y'm0[t,x,|afc|exp(iek)] (14) 

where Ck now are complex numbers with fixed absolute value equal to \ak\ but with 
random phase. One can write more relations of this type. 

Eqs.( [T3| , [T4|) give rise to the possibility of numerical simulation of the processes 
occurring in the Universe after inflation. To do this it is sufficient to take an initial 
field configuration in the form of a homogeneous oscillations plus appropriate random 
perturbations, and evolve this configuration in time according to the field equation. 
The mean value of in this classical system is the same as that of the quantum one. 
It is trivial to derive the same result in theories with more than one scalar fields. 

Eqs.(13|,14) are the manifestation of the fact that when the amplitude of quantum 



fiuctuations is large (roughly speaking, when the occupation numbers are large), the 
system behaves like a classical one. Our derivation gives this statement a precise 
meaning. 



It is also interesting to mention the relation between eq. ([T3|JI^ and the funda- 
mental phenomenon of decoherence of the Universe. In a situation when, as in our 
case, a closed system starts its evolution from a quantum pure state at t = it will 
remains in a pure state at any time moment. The system that we are considering can 
be described, in the pure state language, by the operators obtained by Bogoliubov 
transformation f^. However, when the resonance modes become large the modes 
strongly interact between themselves and the pure state description becomes very 
complicated. In rescue, another description comes into effect: the system can be 
considered as a mixed state, i.e. that of an ensemble, of states (semiclassical in our 
case) characterized by the solution to the field equation with initial condition (|I2|) 
with Ck playing the role of the indices numbering the quantum states in the ensemble. 
Explicitly, the system is described by the following density matrix, 

P^J2 l0[^>X,Ck])p[Ck](0[t,X, Ck]| 

{Cfe} 

where |0[t,x, Ck]) is the notation for a (pure) semiclassical quantum state associated 
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with the classical field configuration x, Ck] (the choice of |0[)f:, x, Ck]) is not unique). 
Eq. (|13],|14D ensures that the average of in this mixed state is equal to its true 
quantum average in the pure state with proper choice of p[ck]. It can be shown 
that the same is valid for other Green functions as well. So, for calculating physical 
quantities one can consider the system as being in a mixed state, though in fact its 
true state is a pure one. Note that the time period when Ck are large, Ck ^ 1, but 
still sufficiently small so that the modes are linear, is the matching region where both 
pure- and mixed-state descriptions can be used. 

4. So, we have shown that the problem of evolution of the homogeneous oscillating 
scalar field that has a direct connection with preheating in inflationary cosmology 
is equivalent to the classical evolution of the scalar field with some random initial 
conditions. We also show that the parametric resonance phenomenon provides an 
effective mechanism for the generation of decoherence of the Universe. We have seen 
a nice feature of this mechanism that the transition from the pure-state to mixed- 
state description can be traced explicitly. 

The result of this paper may be useful for numerical modeling of the preheating 
process and hopefully will leads to the better understanding of the latter. Preliminary 
results seems to agree qualitatively with that of the the Hartree-Fock approximation: 
the fall of the amplitude of oscillations is saturated at some finite level. Moreover, 
after this saturation the modes with small k seem to play an important role in the 
evolution. The work is still continuing, results will be published elsewhere. 

The author thanks P. Arnold, V. Rubakov, P. Tinyakov and L. Yaffe for valuable 
discussions. 
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Figure 1: Integration contour in Schwinger-Keldysh formalism. 
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Figure 2: Feynman rules for calculating (j){t). 
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Figure 3: The one-loop graph (a) and the representation (b) of its leading contribu- 
tion. 
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Figure 4: A two-loop graph (a) and its reduction to a tree graph (b). 
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Figure 5: Some tree diagrams arising from the classical problei 



13 



